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I . In trodue l iou » 

The 1 arotropic monel, conceived y not y more Lua. L . u 
occavies ago, has piaye< a roie iu numer -.ca ± - jr^cas t lug far bcysn.c 
expectations. n T hen operational forecasting rth the arocrou c 
model began in the early 1 1 0 ' s , it .as expected ti.at the more so- 
phisticated aroclinic types ^oulo soon replace the former; honc.ei, 
di f f icu 1 1 ies wi th these moae Is , together ■/ 1 1 n t »r tl.-er improv e - 
meats in the bare tropic monel, have caused co.:ti iuetl interest in 
the latter. 

One aspect of arc tropic iro.. ;hich has been Laves Ligate ; in 
some octal 1 is the dynamic stability or nona currents, uuta; ly 
by Keo [i] and Eiiasen [2] for non-divergent fio , an for the 
divergent case y lJ.L in— Nielsen [3] ■ 

In the first pa* t of this paper, the sta. Llity properties oi 
various arotropic monels ill be studied by c o liferent ap 
proaches, namely a fini-e-diCference method an. a finite Fourier- 
series method. A variety of oonal vino profiles ill re con- 
sidered, including single and <. on ie ■/ inu-spec*.. maxima. The 
results of the too methods are compared to one another as .el. as 
to those of previous investigators. 

The performance of the numerical precictiou models, .hen 
integrated over periods of time up to several days, eexs , or 
a en more, has been recently of special interest, Thompson [4] 
has developed a heuristic theory of long-period velocity varia- 
tions in barotropic fio . In particular, it .as sho ;n dm ti i^ 
theory that a do. bie-jot maxima in the kOuaL vi Ac profile 



l 



develops in t >o ox Lhreo ay fro., au i ul^x ui.jgle ^ t. ac- 
conditions seem t:o e ci.arecl.eris t ic of certain t^pcs of atmospheric 
currents shea a lock occurs. This theory is bascc upon the in- 
teraction between edeies and the mean flo/. 

itfiin-Nielsen [3] has recently shorn that the splitting of 
a single isotach maximum in the meridional :in*j profile (single 
jet) into tv;o maxima (clou ie ^t) is a dynamical couseq.ence of 
the convergence of meridional moiaentuai transport ry the etuie^, 
provided there is barotropic instability in au initial flo . It 
may be noted also that the time scale of the changes in the mean 
xonal flor is tho same in oth theories, as may :,e expectec, since 
the physical model is basically the same. 

The second part of thi*> paper is devotee to a study of 
second -order effects of the perturbations on mean xonai flo./ 
which has do k le maxima. 
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t\ » c* h> u. C L ,. 

Tie N/..:tic ‘ i ^ r . 'lu c. i. > 

(or equivalent ajctro^i . ^ ) . iay t l ,r > cLCu 1*. t 

dorm ti .e to £ ] : 

^-(\/ 2 fj+v-v(vV+n== 

whi.ru Ip"' is th - -.tica.: r ..i c., ; 



= U /. VV' 



UtC . i v. 



.iv 



and q 2 = Jd i£- 

^ s 2 dp dp • 

I;; equation 6"= — (^~) 1 

stability, f 0 a btaiuarc, '& oe o the Corio 
c\ ( p j i. aa wi i , ip n ica 1 i .j iCL to. . - . a ci i iu ^ t < ^ 



. ) 



;a it.. .. w _ 



■ i* e a - 



l i i U L . ^ 



trio norinontal “ino. ,xar the live 
applied ( d 00 mi-* ) . iarue 1 y 

\ (.<, y, p, O = Op) \ to , O'. 

By appiyuv, che perturbation ; ouu , toe poo . <. ^ i 



•o 



c -u r.i, a 



(2.1) may a 'iLtten 1 *i the fv 



V 2 lf + u^( vY') +■ i*\p - £p) = rff 

_ _ BW 

ic tne I'trcur a tic atu.'Lt t _ tetioo 



9t 

Ore U ( [O 



ax 

is the m_au aoi.al 



dr 

LA 



3t 

j 1 o c l > , a o 



The per tar hat ion stream functa. 



CO.: 



iru 



y'u, y,t) = fey) £ 



ii. > 



a 



S U O _ L L t U L 1 ^ \ /L . J ) jlicu \ * U ) j ( l i . » l .) u i l 



(D c )[^t 4^] +- \_e - ~ t- fc] 4> = o # 



to 



(*•*) 



t /> x 2 7T 

“L 



(a real quantity) a*iu c, the pnase 
elocity, may *e complex, teat is 

C - C r -i- t Ci (.../) 

!■, ia (2.7), Q the. amp .iLuee >i : i e an . /. r uuuut ia i 

.•auction of l uue t an^ the a a is t>ai to e unsta 1c . 

For £; > 0 > Lae ampi •. Ca.e i .. i icreas- ' t h a mu; rf, on 

the other nai»u, C; < 0 > the a v e ill a ampu , 

Equation (2, ) i: a . ecouh-or er eifier: ntLai equation an, 
lo e solved j.or co.“La»n ounnary coo. . i t iOu.. , o uicc ^ x^tur’ a ce^ 
j-i* - oria i CujlX*0-. wS a, e norma- y Guru 6v, to t n*- aorta au s^/uti , enc 
simp lest type of o. La le ai rn>ar\ c JiriL'o u ay a ,riuc i iu 



i,;e t.orm 



(f)( 0 ) = (|)( D) = 0 , 



(z.S) 



cox iu:- pon m 0 Lv. x j oUsK r. r j.l at j — o a a y -uy ..ttut iuu 
per tui : a t ion stream "unction i.*, a^u.nee to vauriu 

No;, ciu an e‘. yn-,a : u^ pi o- .. e-.,, .. are interests* iu the 

reiat iousbip et\ tea c, yn auo ot* er per Liu*;:. it parameters fen 
hich there ext urn. Lions .h tne system (i.t) aac (2.8). 

Since tne a sic err rent is a function of \ , tne mathematical 
pro'. Lem here is rather cornplicateu ; a»r generally re cannoL fini, 
vu for s .Maple ranci ie..o U (y), an analytical expression for the 



pcriiii ation amplitu .e e (>). Hence it appears logical to attempt 
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solution » y a 
series met not 



iriii . u'. * . j, : ru'inct. mcLuo or a tl.ti- - Four .ex 
. iiicu ci; i aua^te- Lj uii ^arl ouiar cun. a 



cone it ions (2.o). 



B. i'lviUliO*. O... ... ..tit:. # 

In the a t.-rcU . [ U^.rc. , xociij, p/_> . ll ~j c . , tci e. 

^cujrait} na\e ab; , cci.a. - .^lr ; u.i c-, r .c. t>c...e lacitu w 

For simplicity, particular syiv etric p'o'> lies, c-sp^cia '» ly siu 0 c- 
^.ei and oou re-jv: proiius, /u. ie ..iscaa^u , 

Since only s y t..:;^Li .c .inf pi or ire'. ane ou.r aiy, con - it ions 
are consi »creu, syi.ruetr;, o . the arrpiitu. e 0( y ) o; t: e perturl a- 
t ion stream faiiCLion ay o u,q rcte To ere . . . u, July iiai ‘ o 
the nine belt .e.u e Cv/siotre 

lu so l oinj, e ibau consiuer she iruncticn p(y) as ^iven 
y its values ^ at In points in cue incur .'a I, 0 ^ y ^ 0 
such that the iu^nv . 1 is i\ •_ eh y these pcjiLo nto N + j. 

uqeal parts, eqch o the. icn i.n ; = u (figure 3 . The 

N+l 

meridional eeri.ati.u^ 01 p(> ) •■ill c approuinaiuv -y the cor- 
responein^ quotients of eenterco * iifereuces seen as 

fi!£) = iiM + 4vi ~ 2 fe 

''dr'” d 2 • ( "' 

In general it may. e e-q;eci:eu tnat the a pproaii. a t ion of he 

..eri vatives ■ y this c* if ference ratio jr-.c ectej. results 

as the value of N increases, althou 0 h the mathematical, analysts 

. ecoaies more tedious. Ho. ever, _or simplicity in this section 

only, the cases of N -- j am. N - 7 .ill e consiverec. 

Tart I. 15 - o, (or =- ). 

o 

By combi iiri,j (2.3) aik. the ji\eu .Quinary coauilious (2.8), 
the aeriv{| lives of the perL^r atiou a. plituuc at each point 
l '2 • i) can he written in the form 



o 



(2.10) 



(£t), 

\ d. 



2 <f>, - 2 to 

d T ~ 



<fc -h <1 


>o- 


■ 2 '1 




C 




h 



mere L no ouL-scr . : .;li 



- s v- 1 a. l i c va .t at 



D _ 2 

2 5 7 " 3 



-J , anc v = Jj, respectively 




Since (2,6) must ’e satisfied at each uoi..it o 



1, the 



.?ollo.;ia 3 set ol equations uith the use of (2.10) are obtained; 

[(o„-o(-z-/j ! ) +di^-D; , ) + .i*qx]t + 2(ti.-c)i P = o 
(iJ,-c)-h + [(n,-C)'-2-H i d z ) + di | S-II,'')+d"fc]f> 1 +(tIrC)!h=0 < 2 -; 



here U , U, , U,. , U , D. 

l'l I J (\ 



auv 



are CiK: value:-, on U a,u U 



at v 5 y ana y' , respecti\ ei\ . Here the i : ou : ;e prime rears - 

" O i /L 



seats the secon 1 derivative with respect to y « For a uno jii 



initial venal \/inf profile, (/. .11) cons 1 1 taiies a system of three 

equatio is in three . :..io'c',is $ , •/>-,> a. a o,. . Si.icc Cue sun 

on 

(2,11) is homogeneous > the posoi* J.e values of c may no - la 
obtained from the condition that the ueteriitinanL of the sys \ 
(2,11) mus t van i s u in ore) e r to i :a v ... no n - 1 r 1 v la 1 c > o 1 u L ions . T his 
conditions leans to the so-called frequency cq ation in the form 



c J [ 4 ( A - 3 )] +- c [* T f < /u Z H *- (v, +- d 2 ) ( i - ^ ) } 

+ (^ -/ ) ( Fi-D. >')+£/? fu ,+U, ) - ': ( li-Hj r)]+c [/n n, Jj z (V 
-/ ) -t-n II -u u, +15 U, ) f +• i A (fi+ft ) +(u, +TI 2 ) f / - r /, ) 1 • 



./here 



r\ 


- h U.7 J|T 


(It j , ; ( 


ii -u..r)] + [; 


n.iiztr 1 /) 


-r( 


l'i dz~*~Tz U; ) r 


i ) -TT,r. - 


2 TJ. D, ( u 


r)j = c 


i o 


= d (3 -Vo 


b , 






?, 


= dhh-Ti ( " 


) , 






?, 


= /(.d-U. 


) , 






r 


— 2 T- O'd ■ 


•> 







and '^-y+d'T 



It should be note that <£p aau 9 .ere aioume constant in 
deriving (2,12), Since the freq.eucj, eq .atio ; (2.12) is a cu 1 ic 
equation in the phase /elocity c, iu ge.^ral tnem e:;ists throe 
roots, which may include a pair of complin roots. The imaginary 
part of the phase , ^iocity, aich i necuojar^ for unsta' ie 

perturbations, .ill ,.e consiaeres of primary interest in this 
paper . 

The frequency equation (2.12) has been ,-ofvc. numerically 
for a number of cases .ith a Control -Data-Covporai-iou 160n 
computer, using tue so-cailee Muller metnoa for finning the 

roots of polynomials. For each of the initial /im profiles 
an:, meridional f hLs cancv-s i ctvec.i the ri^lu ouudaries, the solu- 



tion ..-as obtained. for various renal .a' ele :gths (L 



2n 



JA 



), 



y 



t-.ee pm^ ^ ai. o cj 11 . Ther tie ect . o m y, a»iu 
/ / 

divergence q ere teste iy making each o f thtuo parar.K tors :ero 
in a series of computations. 

It is of special interest to consider the stability problem 
for those types of initial wind profiles which arc not onLy sym- 
metric with respect to the center of the channel, : ut also ire - 
quently found in the actual atmosphere. 

Example !• Single -jet asic current. 

For the first example, the . ona^. wine profile /ill be con- 
sidered in the form 



U= B (J -cos 



2 ~ 

V 






(2.13) 



which is shorn as curve (1) in figure 2. Inserting (2.ij) into 
(2,12) and varying the pertinent parameters ;3, E, , a. 

q leads to the results given in figures 3, , a, 6, ana 7. 

In figure 3 the meridional width of the channel ; D, is 
2000 km and the maximum uonai speed at the center of the channel 
is 60 m sec ^ (E = 30 m sec ^). This case corresponds to a 
rather intense single-jet /ind profile, ..nich is often found at 
the 500-mb level. The a scissa in the figure is ^onal ^avelengtn 
in units of thousands of km, while the ordinate is the imaginary 
part of the phase velocity c^ . The figure contains five curves 
corresponding to the cases: curve (i)^"- 0, q = 0, (2) p — , 

q = 0; (3) p- pt,? , q 2 = 0.5 x 10 1 " m “ ; (4) (3 = 

2 - 12-2 2 -12 -2 
q = 1.0 x 10 m ; ane (3) 5 — , q = 1.3 x 10 m 

Here represents the meridional erivative of the Coriolis 

parameter at Latitv.ee 41N which has the value 1.618x10 ^sec . 






By comparing L Irst cur . , _t .oc. t at ihu para- 

meter £? has a ataM . influence u i t e r rt .r atiun jl the 
sciiSe that Li.u ma /..icu' u -f c. is smaibe.. a t .e ai. of ,a\,- 

longths for ./hie. i.r- ta..ility occurs is sure lasiO for j3 ^ 0 
than in the case of ^ = o. horeover, .// conpari .g the first tu 

Curves v/Lth the last three carves It is aLo seen that q, which 

is a measure of the intensity of the divergence Liiplied hy this 
model, has further sta irining effects. Tuerefocu, ine divergent, 
oue-parameter model is more sta le tha the nuii-di T urgent model. 
Furthermore, the stabilising effects of and tne divergence term 
are much, more pronounced for relatively ton or /ave lengths in the 
previously unsta’ 1 e ham's. Thus the maximum c_. ' s occur at shorter 

wavelengths for jS ^ 0 and q+0 , than in ti e case of ft - q = o. 

It should w e note' i., this figure tnat the instability occurs 
only for the intermediate ..onal ? avelengtks (in this case 2000-0000 
for ft = q = o), which is in agreement ith the average length of 
observed synoptic waves. In other ;ords, for sufficiently short 
and long wavelengths, instability does not cxisc. This is similar 
to the results of the two-level baroclinic instability discussed 
y Thompson [?] . 

The preceding results are nicely verif ied h ; fig. res 4, 5, 

and 6, which contain similar computations f it i thu meridiuiai 

widths, 3000, 4000, ar.u 6000 km, respecti /eiy . The additional 

information shea, in these figures is that the banc, of unstable 

wave lengths is shifter to .*aru longer .a\elcngti*c such that c^ 

becomes maximum at £ 2.0 for & - q = o ; an also the 

D 
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ma^aitune o t i v‘ ?• c v t c.-a- u i ; i La,., ca^ . 

TniU result i.b corvla :. . c.. i t it n o 1 : .c : [ r ^ 



it may .o j.Ol(\. : aim the ;a ,11 i,.^ ±ecto o a a*- q arc 

vrerally more proaoiiiicc' ? itn increasing D; t .vs the a^es arc 

more stable in the hireryent model - it.i v icrcacV. ; /idth or the 

channel. For example, *ith u - 000 vm, al. : :r th^tiis arc o^a r.' 

in the case oh 3 — pW a<1 4 = 1*5 - 1 0 \ anc only the 

arcs in a narre ■ -ane arcc.r- 5000 v.r are emtu ir for ft — 3^ 

2 -12 -2 

and q = 1.0 x 10 . . For b - 6000 in, x ry the case, 

6 = 0 and c. = 0, io u.ista: ie. 

Figure 7 contains five curves ^ivin 0 the ima jinar^ root c^ 
as a function of he meric ioual /icth ? u, "or a fixed value of 
the aonai v/avelenyt , L ■- 5000 kin. 



sample 2. 



For the secon. profile, it oili he assume- inat 

Z = 2B 



u i - 



u = 

c- 



( 2 . 1 :) 



T itcrc B = 30 n sec 

Usiny the f initc-c.if f erence approximation for the derivatives 
o.ial niucl spec ^ :ith respect to \ yielos 

rr" - 

(Jo = n — 

o 

rj"— _2L 

U * “ 2 cP 

U ,"= 0 



or 



( 2 . 12 ) 



li 



Inserting 3,-u* (2,1.) l j ' K . . ) ci ^-r.o,. i.ig lI.c aucc , '.ar 

computations leaf s to ti.e res ets ho..i» -' i;ur-: 3 t, 3, an 10, 

corresponding to L 2000, nOQC, and 6000 km, respectively . 

It is readily seen that the results of the instability 
feat res for this profile are q.ite similar to those of the first 
profile, which is to e expected since the to profiles are similar 
By comparing with example 1 , the main difference in this -inn 

T 

profile is that the nan imam 'a lues of c. occur at — - r.J, foe 

1 u 

D > = q = o, and also the ana of unstable ./a\ e lung ths shifts to 
slightly shorter .avelcngtas than for the previous >3sic current. 
Moreover, when b = 4000 km, it is seen that all aves are stable 
iii the last two cases. Hence it may be coacluocd that the sta il~ 
ixing effects of and q in this basic current are much more 

pronounced with increasing D than in example 1 
Example 3. Dou '.e-jet asie current. 

In preceding t o examples, a stability stuiy as given for 
syrmuetric single- jet -Jim profiles* Now in this example, ..e wlIjl 
investigate tne sea* ililjr character Lies of the symmetric ^ouj!.c- 
jet wind profile, represented ■. y the expression 



U-6 (si n -yV ) (2.io) 

ahich is shown as c rve (2) < 1 figure 2. 

Inserting (2. Milo (2*1^.) ana perform! .> rue series o" 
computations yielus the results shown in figures ii and 12, cor- 
responding to D = 20C0 ana 4000 urn, respect! . el,. . 

Contrary to the results of example 1, it is seen that 
ano q go not show the stabilizing influence for all wavelengths. 
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F^: relatively ioha:*. a.ciui/u.b (ic.., „r l» an uit a . elen ;l,. at 



..hich the maximum of c. ..or p ~ q 


ss 0 


OCCUi J ) , 


the 


^ — (2 fleer 


'nas a uestabiliJ Lug influence on the 


per 


turba Lion 


s i.l 


s xh a ..ay 


that the magnitude of c. is greater 


for 


& 4 = 0 


than 


for =5 o 



Also the i ana cf wavelengths for which instability occurs is 

reader for j3 =%=. Q ■ Furthermore, for these longer wavelengths 

the effect of q, a measure of diverger.ee, is to further increase 

the instability of the . a* es similarly to the ^—effect, 

/ 

It should be noted that tv/o maxima of c. occurs Cor 
D ss 4000 I;m, ./hen f 8 0 am' q rjo o ; and that the waves are 
much more unstable -hen the meridional . ic.th jj is increased. 
Hence the divergent, one-parameter mouei in the uouhle-jel wind 
profile appears more wwista' lc than the non-di .xr gent model, .hich 
is opposite to tea results for the single-jet profile. 

It should be mentioned here that the soldi. is ion of the 
basic channel finite differences into 6 nosh lengths cor- 
responding to the case N = 5, is proablj a cruder representation 
for the dou'le-jet asic current than for the .ingle-jet- hence 
the results for the former may be expected to i e Less accurate 
die to larger truncation errors. This aspect ill e considered 



again in Part II. 
Example 4. 



U=B(M.S^-C0S^g) . 



P (! D 

This profile nao easterly wind components vO tue nortn and 
south of a belt of westerlies in the center of the channel. 



(2.17) 



(figure 33). The values of c^ as a function of the zonal wave- 
length in the servo; of computations are sho :n in figures 13, 



L. 



]>:• , and j.j 



C r.£ ^ 



Cuu I.. 7 re. 3 i>ccLi. .el y 



c ^ y — 000, f 0 0 1 j e 

It se*uu that J am* q ha* c a jta' il‘ ' u 0 influence as 

in example I, ut a»* uasta le a a-.; of . a-’nicu^ths -oes exist for 

the case :?= ^ a. a! q% 1 . l 0 n ~ , even for D = 6000 kn. 

The maximum values of: (for yS ~ q - o) occur nearly at 

L = 1.0 anO are ^rcater than for any of the profiles previously 

»’ iecussec * Ho* ever , the Uau of av cleave he iu the case o 

j3 - q = o, fen. .hicl; iastahilit) occurs, L: :\ud .*arro .*cr tla » 

th-o others „ 



It may e co.xlu cm t..at this profile i.. the ^iverpent 
he more unsta-Le t. an the oir.^le- jec profile in example 1. 
rart II* K - 7, (^~~ ~ £ )• 




mo 



1 



The general eall iir of thu cone asion.; hich have been 



readied in precenin*. parts nay e q.estiona* ccu.se o i: the 
finite -difference mesh-si :e co responding to h - 5 . 

As a check, on the previous results, it ill 1 e of interest to 
study the case, U ~ 7, as epictec in fl.yure It* In this case 
the frequency equation (2*12) is replace. \ a "our th-degree 



In 



equation of the f co»i 



C^h/h ) T\j + C*\ {-Pi £>,+(* o$>o) 4- /?V 0 j + iV (Ao„ 4-Go/Si 
4-2V, J 4- „ jVlS=? r J7 0 S? — 2V? j == 0 # 

The notations in (2,18) have the following -.ef i.nit ions : 

?r -h(^-U ; ) , 

he = f. -U P T , ( // = ;,/,?- V • ) 

| = 2 ■*“ // d 7 
/< - ) + d ^ , 

T 0 =f"-I , 

17 = AtC- 3 H~) -H ■d 2 +-u 3 ) ; 

s 0 =/-n>/) 

S, - A (T, ^Ui 4-1V ) +- i,T -$.< , 

=■ [ (11 — TJ, u_ ) +6 ,n 4- ( U, t U- ) , 

s, -fr.T* - r/ 0 ,u , 



(2.13) 



V4 - T- (T +-1, ) - T, , 

V * — T| i •> 4- L’i ) — mil i o — Ta •> 

V* = 71 (lI 0 ^Ui V-IUYT, . a. id 

V 3 ~ l *& T3 1 j ? 



1 j 



Solutions o ati.o.. (2, jlo) uavt .'cun o. l a i .. l u y the Nullcr 
rr/jf >d [6] for tne sa.ac: scries ol examples as for tre case, h = s . 

Since there exist differences in this., iu.ta i. ity char- 
acteristics of single- jet ana tou le-jet profiles, it is of 
special interest to examine these cases usin^ the finer mesh- 
sice corresponding to N = 7 . 



Example I. , Single- jet basic current (fi . 2, curve (1) ), 

U = b (I -COS^-y) (5a a., (2.13) ). 

The values of c. as a function oi ..avelenth in the various 

i 

series of computations are shown in figures 17, 18, and 19, 
corresponding to J = 2000, 4000, and S000 lea , respectively. A 
comparison to figures 3, ^ , ana 6 sho /s the two cases, N = j 
and N =7, to 'e -.cry similar. However, some differences do 



exist for N *-= 7 as follows: 

a) p has a f.es LaHlicing influence for relatively shorter 
wavelengths ; 

b) the band of unstable wavelengths is slightly wider; 

c) the maximum values of c^ for j3 = q = o occurs at 

L =2.3 iuGtcacw of at 2.0; 

D 

d) due to the lesser sta! Hieing effects of $ and q with 
increasing meridional ./icth, D, the unsta le curve (9) 
in figure 18 for D = vOOO bra cine curves (?) an*.; (3) 

for D = 8000 urn have appeared (figure 19). 

Example 2. Double-jet asic current (fig. 2, c.r/e (2) ), 

U = b (Si fi | V V * i t! ^ l J ) (Sa.»c as (2 . 16 ) ) . 
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A.i? inoicat- : a .j. j. e r , i e a t i n jl f f e r e j i g ^ s ' u jl . ^ u t v o 1. 1 nave c e . 

c.»wCCtGu in tnis case i Lcr,‘ioG cue (J rtatoi* coi.ip lenity o the sou io~ 
jet profile is not as -ell i eproiseiiten l y the subdivision into C 
parts (N = 5), 

The reSw. Lto r or M = 7 is this profile, she in figures 



20 


and 


21, are also q< ite similar to those for J 


- 5 (figure 


11 


and 


12) except as follows ; 






i) 


when D - 1009 4 ; ; 






a) 


the magnitude o 1 c . is nearly double. : 

r 






b) 


since the ■, esta: Hieing influences of 3 


and q are 






much greater, tie ana of uasta le a*, ale. 


i^ths nas 






wieeneci ; 






2) 


w ien D = 4000 la , 





a) tie maxima^. value oC for jS = q = o occurs at 



L - l*o; 

D 



tire last 


three car - ' es 


(3, 


4 ana 


of fi. 


; urc 21 no*, s 


ho . 


a staMl 


iw:.i'5 effecl . 


. ith 


ine real -* 


* , q 


hr wavelen" 


the 


jiven in 


this r i,ure } 


h'Ut 


i l bC ems 


*P o' o a 


le tna t the 





effect of q /ill he reverses for sufficiently longer 
..ave lea j ths ; 

*■ 12 - y 

c) for a.- q fc a l.j 10 n ^ , c. vanishes 

at L = ^000 ha. 

In summary , it ma) i. e stated that the essential features 
of the stability characteristics are similar for lue cases N = h 
anu 7 ; however some details appear to c lost in the cruder 
sJ. division of the asic channel corresponding co the case, N = o. 
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C, Four ier-sa/ L . iUuCk , 

If a series of trigorc nctric terns ir use co represent 
U (y) it appears appropriate to seeK a soiuti . A: (2.0) for 
p (y) as a fiuite Fourier series of the form 



M 

= T sm 

n=i 



, M= / , 



/- 

J > -> 



w 



C2.rn 



mere 



X = 



7T 



v • 

Then the expression (2.1J) Deviously satisfies the houndary condi- 
tions at the rigiu ails at y - o and at y s D, 

As the value of U in (2, Id) increases, more reiia le results 
regarding the star iliLy features may 7 .e expected, . diin-Nielson [3] 
has found the stability characteristics for aiues of n up to 3. 

Larger values of h (up to n = 8) have een considered by 
Eliasen [2] but mostly in cases where the effects of and 
divergence terms mere excluded by assumption. U s results sho^ec 
the existence of only one unstable mo^e correspon i ug to one 
values of n. an^ also that c. at t = 2.0 in the same for 3 - A - b 

L) 

but some /hat smaller for u = ? ana 8. Hence thin study /ill c 
concerned only nith snaU aria ooc values o f n. 

To verify further the results in the foregoing section, the 
symmetric, harmonic uonal velocity profiles uil. 1 again he considered 
L), Sin^le-^et asic current: 

(same as (...13) ) (2.20) 

./here B = 30 m sec. 

Considering the case h ~ , inserting ( .20) and (2.1>) into 

(2.b) yields a system of linear homogeneous equations 



U^3(i-CiS ~ 
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[( <n c f P - f ( 3 .tf-X )] p F t J ( UA^-/r /> p =0 

[f r^-3X" ^ +[ff2/xV// i )]<b=<? (2.2D 






la order to have uon-trivial solutions for 0 , 0 , ana o K , 

to -> 

the determinant must vanish, ./hich ^ives the frequency eq.ation 

C 3 ( RW E ) - ctEiRV+S W) (-KW«] + C [t fiV-TH ) f & (K\'l S Ia)) 

-FHtv] f [ Ct ( SV-TH ) - FHV] = 0 



( 2 . 22 ) 



/acre 



E - A *+u z +r , Er = jS - ) , 

H = ) , p = | -(JU Z -5K) ; 

R =- pAvaU4-<r , £ = 3-B( Y.a~+ /,P ) ? 



T 

and V 



“(2 / AV .// 2 ) , W = 2 k A^-k-.U 2 +- <P 
?-&(2(rA-+M') . 



In the case, M = 7, the frequency equation (2.22) is replaced 
by the four th-de^ree equation: 

cnYRWF)tC 3 (YV J fZKWE)+-c 2 (YK,fZK 3 -TRFx;^c[YK, 

(2.23) 

+YK| -(££> HRtr )TX] t [HK 2 -(<tS - pH )TX i =0 
..here X = -f-( M~ k - */-£ A s ) , Y = ^9 A" + ll~ k- 9 i , 

2 = ^-P(f9A 2 a/d) , rio = RV t-SH 7 
M,= SV-TH > K d - En o +RW^r > 

K, =Etl, h-CtHo - pH IU 7 
and K2 =<t‘ v >.-?HV . 
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The roots of (2.22) and (2.23) lor the prt ion series of 
computations are sho n i lt figures 22, 23, 2n , 2^, 2u, aiu 2 7, 
correspouu iu r ' i.j >. ^ 2030, oOOO, am wOUO .i ."o .. the case:;, 

N = 3, a uu N - 

There is o viously a close similar . ty la the essential 
features of tie stability characteristics anouj, the cases, N = 3 
( ,iven in hiin-hielsun [-j) , d = ^ a a ii - 7. Furthermore, 
tiler e also exist s imi iar idles etwcon thu res. it;, of the finite 
difference raethov and those of the Fourier -scries lutliou. Ct 



c t'iOtee that the results of c. e cas 



(N =3), in 



the fonaer metixx are particularly close to the res. Its of the 
case N = 3 in the Fourier-series uutho'. , waile the results or 
the case, a = hr (M = 7), correspond more ueai 1} to those 



It may 1 e concluded that the results for small values of h 
generally represent the essential fuat.ires of the sta ility 
characteristics In the quasi-!. arotropic flow; ho. ever the folio 
iny vtifierences iuh iucieasinj N o appear: 

a) the max .•'■.nun value of c. "or c? = q - o is less 

i ‘ 

(c, - 3.3 7 m sec for N = o a no 3 . u 0 m sec 

i max, 

for N = 7 ) ; 

) tine sta : i*.i. in, crfects of a.i 
are less t . ronounceu ; 



q . jl r n nc r ea s tn 0 u 



c) the maximum \alucs of c^ are shifted to arc ion[;er 
wavelengths sich that for the case ft ~ o ana q = o 
they occur at ~ = 2.0, for H = 3 ; - = 2 , _> , for N = 5; 
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I 



cl ( i U ~ j » C J u i i - . / , 

J 

1 ) the aa'. o; . al /a*e ) eu'-tlii or uica lusta L11 l/ 

occur:; eco; ic . , 

2 ) . Dou: i e* ** j e L ds ic current ; 

U - B( f-C0S&-y) (2.-,) 

;here H = 30 o- .>co 

Thi i profile ha* J = o at otn ri 0 iu alls and at the center 
o r the channel, and also nas u - 2B at y = t an at y = t h 
(fin* 34, U (o) ). 

By inserting ; a.r , 'A's) (1\ = . ) into (2.c), the folio *• 

in,-; set of linear : roino 0 ou-» ■ .t. : , .atiods Is o’ taL* • : 



((.AV/'-fCHc f ;3-3(A 2 +A 2 '] r, + [f('lK-jr)]^ + ( f- (JK+/I 2 )] 0 

[f(/O0-//v]'i>, t((9AVA'-'fq i .)C f p-PO^+H* ) } = C ( 2 . 2 .) 

(f Jr A* )] P, 4- [ ( Ztf+M^q* ) C + 8 -B (ZtA* 4/^)J'hr = 0 . 

In oroer to have ao.:-:;ero j-o*. /^, sne the 

determinant must anish* ih L con hio i , r us the frequency 



equation 

Cqwt'RJt-CqWiE^tR'j) n'Eh . + 1 IjTR tW . n -pH) 



A\£ 


S + KC-'j +[TTS+V('6>:- 


_p. 


j 

/ 

1! 

o 




here 


F = 






= f - f : 


+vn , 




H = 


F /■» \ //> > 

— i / • ' M / > 


T 


- 9 k \ 


// - ) , 




• — 


| i , 


h 


= 1 /Vi f<~ 


-+ j ~ , 




s — 


p K ^ 9 A +- - i 


to 


= aFa"-(- /■ 


•> 

"-or, 


and 


V - 


3 -p ( V : -f ) , 









(2.26) 
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Tbit; freq..e;e;, cCj.ati.OT. (2. .u) ,ia;> u .1 .tu-, .ia.> ricaxiy 

for the previous series of compute l io.u. . The values of c. /ith 

£ 

respect to wavelength are sho m in figures 2.o, 2j and 30. 

It is seen that ft and q effects o not 0 ive a stabilizing 
influence for all wavelengths. For relatb eb longer ave- 
iengths, fi and q have a destabilising influence on the pertur 
bations in the sense that the hand of unsta re nave lengths is 
much broader than in the case f3 ~ o or q - a . Moreover, the 
waves are more unstable with increasing meridional iuth D . 

Hence the divergent, one-parameter model in this double- jet wind 
profile is more unstable than the non-divergent model. These 
results are consistent ./ith those in the previous example for a 
double-jet profile in which the method of finite differences has 
been used. 

For this profile, if e v ea values of u (n =2, , 6 ) in 

(2.19) are included, there exists other amplifying modes aifferen 
from those previously discussed, /hich are not shown in this pape 
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3. 



Stability propert s 



fit: stiatifie 



The vorticity equation foi this model has e._., o.-tai.ice l> 
'Araason [8] and may he riti.en in the form 



(v-r’)ff + v-v(vY+n -i 2 v - (3 .d 

w,ci-e tj' = , V = , KxV’/ . “ ad V=Ui . 

Here H is the thickness of the stratified la^e-, a.i i c . ' is tne 
acceleration of gravity corrected for compress i iJLity. 

By using the pertur: ation method, (3.1) ca e reduced to the 

form 



(vM^fT+u^cvVVff^ 



ill. q 2 £J 1 _ 

a !~ 



(3.2) 



here 7 is the perturbation stream taction a ne U = U (y) is 
the mean zonal velocity. 

Substituting 'Vf' / ( t - 7 ■> /* ) — <J> ( 7 ) .£ ^ ^ ^ ^ 

into (3,2) yiel s 

(D-c)fi"-,r<i)-r i 3-o'dq i (c-ii )] h — o (j.3) 

where the primes -e,.ote eif ferentiat ion ;ith respect to y. 

In comparison ith (2.6), equation (3.3) nao one more term 
which is peculiar to thi; moael. Thus it may e expected that the 
stability characteristics in tnls monel ill uiffer some. /hat from 
those of the previous divergent , vue-parare. Lew model. 

The analysis will ! e restricted to tie initial ; oral in, 
profile of the type 

U — • B C I — CO" (1 ) (same as (2.13) ), 
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iL cUjO L j id U 



JL jl GC| Ul.! c y 



■» / v . j ) . The rci .it- 



C-- ( n F ) 4 ' * 0 ( R + S 10 + K U r] 4 c [EC 5. V - ; 

■KriKV+i 0 ) - p|-| '•:] i j & ( $ V-TH ) - j H V ] = 0 



E = X ? +/!= + q 2 
H - -f-(^Ar+- tf-‘+ V ) , 



G = ? - f-(v/V 3 r-v} , 
"K= u' i q'-s v) - 



R = 9A* +/j f* + q* , S = f?-£,;r+,^2h , 

T = f( Z I V 4 M - 1 q ' ) J U) = 2 Kv 4- // " -H q 2 
nd V = -• - h (Zy v + >r ^ r l i 

* 



The res, its of e e t rcujus cun^rtal l • - a sirits arc give.- in f i u ures 
31 am 32, corresponding, to n - rOOO an*. oOOQ :,rru 

It is nec.i y comparing tie cur oss (2) a^d (3) m figure 31 
that q has a aesta. iiizu.g influence f 'or reiati. ely louder 
areiengths (lou 0 c*r than the navels. io tn at •hich the maximum 
for curve (2) occurs). Moreover, tuese effects are erified rith 
increasing \aiues o,. q, as riv^n y curveo (4) a, e Hence 

the divergent, stratifies ii.oi e i is sore unstable char the non- 
aiveryent one. Also comparison of figures 31 anu 22 sho .s that 
the stratified w •..! *s i^.re unsta. -e inn i the iou one- 

parameter monel „ 

These resa-h, arc also veri.ie^ oy figure 32, hicn contains 
similar results :cr ^ = 4000 ku . Note that t e .axes are more 
sta' le wits, increasing u. never the less , lor the case p — > 

q*~ = 1,3 n 10 ~ and u = 4000 l,rn, it seems possi> i.e to have 



d V JL . U 



iti.CH £1 1' C li^t. S ! I • j i l 



cl » iiilu O L 111 is Lei . cUvviS Ci u x. J 

in this figure. 



/c. ~ 



I'j 



i-r . rlome.itUM CraaJCS . 

It is of special interuot o sll nal may v_ inf urn.- from 
the solutions o Lais particular proLie a ai-out l.lu effects of the 
unstable pertur atieus on the basic flo. through momentum transfer. 
Cnanges in the >o.:a current -./ou’c- alter me profile so rauically 
that the initial solution of tne perturbation problem no longer 
applies. This implies tnat the asic profile may, in turn, alter 
the structure of tne perturbations, 

Recently i tne paper y .viin-hielsen [3] , it ,.as eon fou*K. 
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l o .-.e-jet pro'ile in ... »e or t./o nays o. e to mori^ional momentum 
transport \ y the a- fins. 

It is obviously of interest to determine hon an initial double- 
jet min.. profile ill Ci^a as a res. it of momentum transport by 
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From the equation o Motion, the cnauges im the mean xonal current 
can e expressed in tne form 
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Fig. 3. The imaginary part 7 c., of true solution to equation 

(2.12) as a function L of wavelength for U=B( l-cos( 27 Ty/D)) , 
cl = D/6, D = 2 000 km, and B - 30 m sec -1 . Curve (1): 

$ = o, q = 0; (?) : ft - > q = 0; (3): (3 = , 

q 2 = 0.5 x 10" 12 ." 2 ; (4): a = $<t-5 , q* = 1 . 0 x lCf'm' - ; 

(5): = ptf , q 2 = 1.3 x 10' /2 m' 2 . 
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Fig. 4. The imaginary part, c,, of" the solution to equation (2.12) 
as a function of wave length for the same parameters as in fig. 3 
except D = 3000 km. The curves as in fig. 3. 
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Fig. 5. The mag inary part, c., of the solution to eq ation (2.12) 
as a function of v.avelength for the sane parameters as in fig. 3 
except D = 4000 km. The curves as in fig. 3. 
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Fig. 6. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for the same parameters as in fig. 3 
except B = 6000 km. The curve (1) as in fig. 3. 
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Fig* 7* The imaginary part, c . , of the solution to equation (2.12) 
as a function of the meridional width B for L - 5000 km and the 
same parameters as in fig. 3. The carves as in fig, 3. 
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Fig. 8. The imaginary part, c., of the solution to equation 
(2.12) as a function of wavelength 'or U„ - 2B, TT , = B, 

TT 2 = 0.5B, d = D/6, D = 2000 km, and B = 30 m sec' 1 . 

The curves (1), (2), (3), (4), and (5) as in fig. 3. 
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Fig. 9. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for the same parameters as in fig. 8 
except D = 4000 1cm. The curves as in fig. 3. 



• 40 



( i 23TUL ) n 




Fig. 10. The imaginary part, c., of the solution to eq ation (2.12) 
as a function of wavelength for the same parameters as in fig. 8 
except D = 6000 km. The curve (1) as in fig. 3. 



41 



0 1 2 3 4 L (/0 3 km) 

Fig. 11. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for U = B(sin( 71 y/D) + sin(3ny/D)) , 
d = D/6, D = 2000 km, and B = 40 m sec - ’ . The curves as in fig. 3. 
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Fig. 12. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for the same parameters as in fig. 11 
except D = 4000 km. The curves as in fig. 3. 
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Fig. 13. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for ' TT = B(cos . (4jty/D) — cos (27ry/D ) ) , 
d = D/6, D = 2000 km, and B = 30 m sec’ 1 . The curves as in fig. 3. 
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Fig. 14. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for the same parameters as in fig. 13 
except D = 4000 Km. The curves as in fig. 3. 
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Fig, 15. The imaginary part, c., of the solution to equation (2.12) 
as a function of wavelength for the same parameters as in fig. 13 
except D = 6000 km. The curves as in fig. 3, 
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Fig. 17. The imaginary part, c., of the solution to equation (2.18) 
as a function of wave length for" TT = B( 1-cos ( 27 ty/D ) ) , cl = D/8, 

D = 2000 km, and 15 = 30 m sec -1 . The curves as in fig. 3. 
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Fig. 13. The imaginary part, c., of the sol tion to equation (2.18) 
as a function of wavelength for the same parameters as in fig. 17 
except D = 4000 km. The curves as in fig. 3. 
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Fig. 19. The imaginary part, c., of the solution to equation (2.18) 
as a function of wave length for the same parameters as in fig. 17 
except D = 6000 km. The curves as in fig. 3. 
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Fig, 20. The imaginary part, c., of the solution to equation (2,18) as a function 
of v/avelength for U = B(sin( 7 ty^D) + sin(3jry/D ) ) , <i = D/8, D = 2000 km, and 
B = 40 in sec -1 . The curves as in fig, 3. 
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Fig. 21. The imaginary part, c., of the solution to equation (2.18) as a function 
of wavelength for the same parameters as in fig. 20 except D = 4000 km. The 
curves as in fig. 3. 
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Fig. 22. The imaginary part, c., of the solution to equation (2.22) as a function 
of wavelength for U = B( l-cos( \ny /Y ) ) ) , N = 5, D = 2000 km, and B = 30 m see' 1 . 
The curves as in fig. 3. 
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Fig. 23. The imaginary part, c. , of the solution to equation (2.22) as a function 
of wavelength for the same parameters as in fig. 22 except D = 4000 km. The 
curves as in fig. 3. 
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Fig. 24. The imaginary part, c., of the solution to equation (2.22) 
as a function of wavelength for the same parameters as in fig. 22 
except D = 6000 km. The curves (1) and (2) as in fig. 3. 
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Fig. 25. The imaginary part, c., of the solution to equation (2.23) as a function 
of wavelength for* the same parameters as in fig. 22 except N = 7. The curves 
as in fig. 3. 
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Fig. 26. The imaginary part, c., of the solution to equation (2.23) 
as a function of v/avelength for the same parameters as in fig. 23 
except N = 7, The curves as in fig. 3. 
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Fig. 27. The imaginary part, c., of the solution to equation (2.23) 
as a function of wavelength for the same parameters as in fig. 24 
except N = 7. The curves as in fig. 3. 
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Fig. 28. The imaginary part, c., of the solution to equation (2.26) 
as a function of wavelength for ' U = B( 1-cos (4xy/D) ) , N = 5, 

D =r 2000 km, ana B - 30 m sec” 1 . The curves as in fig. 3. 
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Fig. 29. The imaginary part, c., of the solution to equation (2,26) as a function 
of vzavelength for the same parameters as in fig. 28 except I) - 4000 km. The 
curves as in fig. 3. 
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Fig. 30. The imaginary part, c., of the solution to equation (2.26) 
as a function of v;avelength for the same parameters as in fig. 28 
except I) = 6000 km. The curves as in fig. 3. 
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Fig. 31. The imaginary part, c., of the solution to equation (3.4) as a function 
of wavelength for the stratified model, T7 = B( 1-cos (2ny/D )) , N = 5, D = 2000 1cm, 
and B = 30 m sec' 1 . The curves as in fig. 3, 
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Fig. 32. The imaginary part, c., of the solution to equation (3,4) 
as a function of wavelength for the same parameters as in fig. 31 
except D = 4000 km. The curves as in fig. 3. 
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Fig. 33. The function M(y) as defined by equation (4.3). 
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Fig. 34. The profiles of the zonal wind initially and at 1, 2, and 
2 1/4 days computed from equation (4.9) with p = q 2 = 0, L = 2000 km, 
D = 3000 km, B = 30 m sec' 1 , A. corresponds to 5 m sec -1 
(1.39 x 10 6 m 2 sec -1 ), and = 0. TJ (0) = B( l-cos(47ty/D) ) . 
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Fig. 35. The profiles of the zonal wind initially and at 2 and 2 1/4 
days compute^ from equation (4.9) with parameters as in fig. 34 except 
corresponding to 5 m sec“* (A^ = A^). 
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Fig. 36. The profiles of the zonal v:ind initially and at 3, 4, and 
5 days computed from equation (4.9) with parameters as in fig. 35 
except D = 6000 km and L - 4000 km. 
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